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BYPASS TRANSITIONAL FLOW CALCULATIONS USING 
A NAVIER-STOKES SOLVER AND TWO-EQUATION MODELS 


William W. Liou* 

Mechanical and Aeronautical Engineering Department 
Western Michigan University, Kalamazoo, Michigan 

Tsan-Hsing Shih*“ 

Center for Modeling of Turbulence and Transition 
ICOMP,NASA Lewis Research Center , Cleveland , Ohio 


Bypass transitional flows over a flat plate were sim- 
ulated using a Navier-Stokes solver and two-equation 
models. A new model for the bypass transition, which 
occurs in cases with high free stream turbulence in- 
tensity (TI), is described. The new transition model 
is developed by including an intermittency correction 
function to an existing two-equation turbulence model. 
The advantages of using Navier-Stokes equations, as 
opposed to boundary-layer equations, in bypass tran- 
sition simulations are also illustrated. The results for 
two test flows over a flat plate with different levels of 
free stream turbulence intensity are reported. Compar- 
isons with the experimental measurements show that 
the new model can capture very well both the onset 
and the length of bypass transition. 

I. Introduction 

Flows in turbomachinery component passages, such as 
compressors and turbines, are often transitional. With 
the high level of free stream turbulence in the turboma- 
chinery passages, the transition in the boundary layer 
from a laminar flow to a turbulent flow is commonly re- 
ferred to as bypass transition. Bypass transition occurs 
under the influence of free stream turbulence, in con- 
trast to the natural transition which has its origion in 
the development of the Tollmien-Schlichting instability 
waves. The extent of the pre-transitional and transi- 
tional flow region on these blades can be very substan- 
tial at certain operating conditions and it is generally 
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for these operating conditions that the present loss pre- 
diction methods fail. Therefore, a more accurate model 
for bypass transitional flow would improve loss predic- 
tions. Designers, equipped with better loss prediction 
schemes, can then fine-tune blade shapes to optimize 
boundary-layer losses for these operating conditions. 

Concerning with the modeling of bypass transitional 
flows, Savill 1 has reviewed an effort to assess existing 
two-equation turbulence models for the prediction of 
the bypass transition. The models were assessed for a 
few designated test cases, including flat plate boundary 
layers with various free stream TI and pressure gra- 
dients. The results submitted by the participants of 
that effort showed that none of the two-equation mod- 
els tested could predict correctly either the onset or the 
length of the bypass transition region. Compared with 
the measurements, the predicted transition onsets were 
mostly too early and the length of transition region too 
short. 

Computationally, predictions of bypass transitional 
flows were generally performed using a boundary-layer 
equation solver. This is a valid approach and is very 
cost effective. However, due to the parabolic nature of 
the boundary-layer equations, bypass transitional flow 
calculations using these equations have to begin at a lo- 
cation downstream of the leading edge of, for instance, 
a flat plate. Often, not all of the computational vari- 
ables at the starting location can be defined based on 
the experimental data and one has to assume the pro- 
file distributions for those variables. For computations 
based on boundary-layer equations, the computational 
results depend significantly upon the starting locations 
and the starting flow profiles. Therefore, it is difficult 
to evaluate transition model in a fair and consistent 
manner based on the solutions of boundary-layer equa- 
tions. To obtain more accurate simulations of bypass 
transitional flows, a Navier-Stokes solver is used in the 
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present study. Moreover, the computational inlet is 
located upstream of the leading edge of the flat plate 
where the computational variables can be defined based 
on experimental data. 

In this paper, we describe a new model for the pre- 
diction of bypass transition. Considering the intermit- 
tent nature of the transitional flow, we incorporate an 
intermittency function orrection to the eddy-viscosity 
predicted by a two-equation turbulence model. For 
the present study, a CMOTT variable-C^ two-equation 
turbulence model 2,3 was used. This is a variable-C M 
version of the Shih and Lumley two equation model 4 . 
Bypass transitional flows over a flat plate were inves- 
tigated. They were designated as T3A and T3B by 
Savill 1 . The flows are selected to challenge models’ ca- 
pabilities to predict the effect of free stream TI (3% for 
T3A and 6 % for T3B) on bypass transitional flows. The 
computational inlet is placed upstream of the leading 
edge of the plate, where there are sufficient measured 
data to clearly define the inlet flow quantities. 

In the following sections, the turbulence models and 
the numerical platform used in this study are described. 
The results of flow calculations are also shown. 

II. Analysis 
Mean Flow Equations 

The flow properties are decomposed into two parts: 
a mean value and a fluctuation with respect to the mean 


value. That is, 

p = p + p" (La) 

= Ui + «' ( 1 . 6 ) 

p = p+p" (M 

r = r + r (i .d) 

E=E + E' (l.e) 


where p,p, T, E, U, denote Reynolds- averaged density 
and pressure, mass-weighted-averaged temperature, to- 
tal energy, and velocity, respectively. It is customary to 
use both the Reynolds average and the mass-weighted 
average in the decomposition process for compressible 
turbulent flows to simplify the final form of the mean 
flow equations. The governing equations for the mean 
flow may be obtained by substitution of flow properties 
in the form of eq.(l) into the Navier-Stokes equations 
followed by a Reynolds average of the equations. The 
mean flow equations become, 

P,t + {pUih= 0 (2.a) 

[pUih + (pUiUj + r {j + pSi^j = 0 ( 2 . 6 ) 

{pE) %t + ( pUiE + pUi +qT,i + r v u j “h = 0 ( 2 - c ) 


where 

r ij = 2(/i + Pt)Sij — j 

fj denotes the mean molecular viscosity and S,j denotes 
the deviatoric part of the mean strain rate tensor, i.e., 

Sij = \(Uij + Uj,i) - \u kM S X] _ ( 3 ) 

The Reynolds stresses are modeled via the turbulent 
eddy viscosity, pt. In all of the models used in this 
study, the turbulent eddy viscosity is determined by the 
turbulent kinetic energy, k , and the dissipation rate, e. 

fit = C M f„pj (4) 

k and e are obtained from the solution of their respec- 
tive model transport equations. is the wall damping 
function for the eddy-viscosity. 

Byp ass Transition k — e Model 

A new transitional flow model is presented in the 
following. The model is developed based on the ex- 
perimental observation that the flows in transition are 
highly intermittent. Locally, the flow is turbulent as 
the turbulent spots convect through. In between the 
passing of the turbulent spots, the flow relaxes to a 
disturbed laminar state. In this case, a measure of the 
fraction of time that the boundary layer is turbulent 
is the intermittency factor. Normally, the intermittent 
nature of the transitional flow was not accounted for 
in the development of turbulence models. Turbulence 
models were in general designed for flows in a fully tur- 
bulent state. To be able to use a turbulence model in 
the intermittent transitional flow region, we incorpo- 
rate a weighting factor, 7 , in the turbulent eddy viscos- 
ity , p t , obtained from a base two-equation turbulence 
model, i.e. 

Ptr = 7 Pt (5) 

where p tr denotes transitional flow eddy- viscosity. The 
7 function, or intermittency correction function, resem- 
bles the physical flow intermittency factor. Much like 
the intermittency factor that determines the fraction 
of time the flow is turbulent in the transition region, 7 
determines the fraction of turbulent eddy-viscosity to 
be used in the transition region. Therefore, 7 should 
also vary monotonically through the transition region. 
Based on the same reasonings, Yang and Shih 5 devel- 
oped a 7 function using the shape factor of the bound- 
ary layer, i.e. 

7 = 1(H) = 7o + (1- 7 o)7T (6) 
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where 


Hi-H 

Hi-H 


Hi denote the value of the shape factor for the Blasius 
laminar boundary layer. H t is chosed to represent the 
shape factor for the turbulent boundary layer. 70 and 
n are model constants. The model was found to work 
well for the bypass transition of fiat plate boundary 
layer. The model uses the shape factor to characterize 
the state of the boundary layer. In a more complex 
geometry, such as a turbomachinery blade passage, the 
edge of the boundary layer that is needed to determine 
the displacement and momentum thickness is not well- 
defined. This sometime undermines the model’s capa- 
bility in transition predictions. 

In the present study, the intermittency correction 
function 7 is defined in terms of the streamwise varia- 
tion of the peak disturbance energy. The disturbance 
kinetic energy increases as the flow evolves from lam- 
inar to transitional state. Its local peak level also in- 
creases. Since the free stream turbulence influences the 
transition in the boundary layer mainly through the dif- 
fusion process, the higher the free stream turbulence, 
the faster the disturbance energy increases. Therefore, 
the disturbance energy level is an appropriate param- 
eter to characterize the flow transitional region. The 
intermittency correction function we proposed in this 
study involves the peak disturbance energy level and 
the local free stream turbulence level. 


7 = 7o + (l-7 0 )/ 

_ kt-kt 

*; t-K 


where 


*n = K{x) = roax(& + (x, y))| x - fc+(x) 
ke(x) 


k?( x ) = 


«?(*) 


( 8 ) 

( 9 ) 


( 10 ) 

( 11 ) 


model. Therefore, the selection of the value of k+ is 
consistent with the base model. is determined by ex- 
amining fully developed turbulent boundary layer and 
its value is set at 4.5. The term k £ —k+ in Eq.(9) is used 
to normalize the term k+ — kf , such that the transition 
progress variable, /, satisfies the equation 0 < / < 1. 
The effect of the local free stream turbulence is ac- 
counted for in Eqs.(9) and (10). In Eq.(10), k+ rep- 
resents the difference between the peak value and the 
free stream value of the k+ profile at a streamwise lo- 
cation, x . Eq.(10) allows the flow transition to proceed 
according to the internal peak level of the profiles . 

As was mentioned earlier, the base two-equation 
turbulence model was that used in Liou et al . 2 and 
Yang et a/. 3 . The model is described briefly in the fol- 
lowing. 

The model equations for k and s are, 


pk,t + pUikj = [(p 4* /it)i i],i - puiUjUij - pe (12. a) 


pe.t + pUi£ t i = [{fi + —)£,;],,■ - Ci^-ptijUjUij 
<T C K 


where 


- C 2 f 2 p— + viitUij k Uij k 


Ci = 1.44, C 2 = 1-92, <t £ = 1.3 

/ 2 = l-0.22exp[-A 2 ]’ R t = - 

0 ue 

1 

C M = 


A 0 + A,U^f 


where 


= JSijSij + ClijQij 

Qij = &ij — €iji z U>k 


( 12 . 6 ) 


(13) 


Qij is the mean rotation rate viewed in a rotating ref- 
erence frame with the angular velocity u;*. The param- 
eter A s is determined by 


k e and ti T denote the local disturbance kinetic energy j 

in the free stream and the local frictional velocity. A s = \/6cos<£, <f> = -arc cos (VEW) 

The model coefficient kf relates to the threshold for Si S kSki - (^) 

the formation of turbulent spots near the wall. Shih W = — - S = y/SijSij 

and Lumley 4 have shown that turbulent energetic ed- 
dies reduce to Kolmogorov eddies near the wall, where The damping function is defined by 
jfc+ = 0.25, and all the wall parameter are characterized 

by Kolmogorov microscales. An estimate for the spot = [1. — exp(— (a\Rk + <23 iJ* + a 5 /Z/c))]^ (15) 

generation threshold level in terms of is thus set at 
0.25. It should be noted that the wall boundary condi- where 
tion for k is also determined by the Kolmogorov behav- 
ior of near-wall turbulence in the base two-equation a\ = 1.7 x 10~ 3 , a$ = 10“ 9 , a 5 = 5 x 10~ 10 
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Rk = 


pyfky 


Note that the value of C„ is bounded by 0.09 in the 
current application. The near-wall boundary condi- 
tions for the turbulent quantities are determined by 
examining the Kolmogorov behavior of near- wall turbu- 
lence proposed by Shih and Lumley 4 . They have shown 
that energetic eddies reduce to “Kolmogorov eddies” 
at a finite distance from the wall and all the wall pa- 
rameters are characterized by Kolmogorov microscales. 
Therefore, an estimate can be obtained for the turbu- 
lent kinetic energy and its dissipation rate at the loca- 
tion where large eddies become Kolmogorov eddies by 
using both direct numerical simulation results and an 
asymptotic analysis of near wall turbulence. According 
to their analysis, this turbulent limit point is located 
at 

ft = f 06 ) 

«T 

At this limit point, 


k„ = 0.25 u? and e n = 0.251-^ (17) 


In practice, the boundary conditions are enforced at the 
wall. With the application of eqns.(16) and (17), the 
turbulent time scale near a wall, similar to the velocity 
and length scales, is determined by the Kolmogorov 
timescale. Therefore, there is no unphysical singularity 
in the current model £ equation. 

The present intermittency correction function tech- 
nique characterizes the transition progress by the vari- 
ation of the disturbance energy. Therefore, the inter- 
mittency correction function does not necessarily cor- 
respond to the physical flow intermittency factor. 

In the following, the numerical solution procedure 
is described. 


III. Numerical Solutions 

The Favre-averaged Navier-Stokes equations and 
the model transport equations have been solved nu- 
merically by using the COMTUR code developed by 
Huang and Coakley 6 . Briefly, it uses a line-by-line 
Gauss-Seidel algorithm and Roe's approximate Rie- 
mann solver. Yee’s MINMOD TVD scheme was applied 
in all the computations. The mean and the turbulence 
equations are solved in a sequential manner. All the 
calculations have been carried out with the same ini- 
tial and boundary conditions, including those for k and 
e* The computational inlet is located upstream of the 
leading edge of the plate. The inlet flow profiles are, 

U = U u V = Q i k = k l% e = €j (18) 


where the values of U\ and k x are given by the mea- 
surement. The value of is obtained by examining 
the turbulence model transport equations for k and e 
at the free stream, i.e., 


„ dk e 

U\ -7— — — Q 
ax 


( 19 ) 


- C ’t 


The value for e x is determined such that, with k x and 
€i as initiatl values, the solution for k e obtained from 
Eq.(19) agrees with the measured data. 

The computational domain was covered with or- 
thogonal surface-fitted meshes, with grid clustering 
near the leading edge and the wall. Several differ- 
ent meshes were used to ensure grid-independence, 
by varying both the number of grid nodes and 
the grid clustering. The results of the grid in- 
dependent study are described in the next section. 
The results obtained by using the Launder and 
Sharma (LS) turbulence model 7 , which is regarded 
as one of the models that can predict well bypass 
transitional flows, are also shown for comparison. 


IV. Results and Discussions 

In this section, the simulation results for T3A (3% 
TI and zero pressure gradient) and T3B (6% TI and 
zeor pressure gradient) are reported. 

• Fully Laminar and Turbulent Boundary-Layer Flows 

Since the transitional boundary-layer flow over a 
flat plate consists of a region of laminar flow before the 
flow become transitional, it is essential that the numer- 
ical solver be able to calculate laminar boundary-layer 
flows. To validate the Navier-Stokes solver, the lami- 
nar boundary layer on a flat plate was first calculated. 
The flow conditions correspond to that of the T3B 
case. In this calculation, the turbulence model equa- 
tions were not solved and the eddy-viscosity was set 
to zero. Therefore, the computational results represent 
the numerical solution of the Navier-Stokes equations 
for a laminar boundary-layer flow over a flat plate. 
The calculated and the theoretical skin friction distri- 
butions are compared in Fig.(l). The good agreement 
between the analytical curve and the computational re- 
sults shows that the solver can faithfully describe fully 
laminar boundary-layer flows. A flat plate turbulent 
boundary layer was then calculated using the LS model. 
For a fully turbulent flow, the model transport equa- 
tions for k and e were solved and the resulting eddy- 
viscosity was used in the mean flow calculations. The 
computed skin friction coefficient and the Van Driest 
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II correlation for skin friction for turbulent boundary- 
layer flows are included in Fig.(l). For a fully devel- 
oped turbulent boundary layer, the results show that 
the LS model slightly under-predicts the skin friction 
coefficient, which has been observed in many other cal- 
culations of such flows. The fully laminar and fully 
turbulent flow calculations results suggest that the nu- 
merical solver is accurate and the LS model is correctly 
implemented. 

• Free Stream Turbulence Decay 

As was mentioned earlier, uniform profiles for the 
mean and the turbulence quantities were used at the 
computational inlet. The levels of the turbulent ki- 
netic energy and the rate of dissipation at the inlet 
were set such that the decay of the turbulent kinetic 
energy, according to Eq.(19), would agree with the 
measured data. It has generally been observed in ex- 
periments that the location of the onset of transition 
moves upstream as the turbulence level (TI) in the 
free stream increases. Therefore, it is important in 
a bypass transition simulation that the computed de- 
cay of free stream turbulence matches with that of the 
measurement. Fig. (2) shows a comparison of the com- 
puted streamwise distribution of the free stream tur- 
bulent kinetic energy with the measurement. The com- 
puted results are the solution of the Reynolds-averaged 
Navier-Stokes equations and the LS model. There is a 
good agreement between the computed and the mea- 
sured distributions. In fact, it is found in the grid- 
independency studies that the free stream turbulent 
kinetic energy distribution is rather insensitive to the 
grid point variation. 

♦ Computational Inlet 

Fig. (3) shows two distributions of calculated and the 
measured shape factors using the same grid of 120 x 120 
and the LS model. The computational inlet was located 
at the leading edge. The distance from the leading edge 
to the first grid point downstream, dx o, are 10“ 3 m and 
2 x 10~ 5 m, respectively. The inlet turbulent kinetic 
energy and the rate of dissipation were determined by 
Eq.(19). The inlet velocity profiles were taken from the 
Blasius solutions, i.e., 

U = U B iasius.V = 0. (20) 

The large changes in the shape factor near the lead- 
ing edge in the case with a small dx o indicates that 
there is an adjustment of the flow from the given inlet 
profiles to the Navier-Stokes solutions. It seems that, 
for the case with a larger dx o, the flow self-adjustment 
has not completed before the predicted onset of transi- 
tion. Therefore, with the computational inlet placed at 


the leading edge, it is possible that the computed flow 
is not entirely free of the error introduced in the nu- 
merical inlet conditions before the transition onset and 
the prediction of transitions can be inappropriately af- 
fected. Consequently, for numerical simulations of the 
bypass transitional flows over a flat plate, it is more 
appropriate to place the computational inlet at the up- 
stream of the leading edge, where the experiment can 
provide data necessary for the determination of numer- 
ical inlet conditions. Fig. (4) shows the comparison of 
the calculated skin friction coefficient with computa- 
tional inlet located at upstream (xo = —0.017m) of and 
at the leading edge (xo = 0m), respectively. The grid 
atop the flat plate is the same for both cases, which is 
120x120, and dxo = 2 x 10" 5 m. The results show that 
there is a significant difference in the predicted location 
of transition onset and transition length. Therefore, it 
is necessary for the computational inlet to be located 
upstream of the leading edge of the flat plate. 

• Grid Sensitivity 

Orthogonal H-type grids are used in the present cal- 
culations. A typical grid of 160x120 is shown in Fig. (5). 

In the cross-stream direction, the value of for 
the first grid point away from the wall was set to be 
less than 0.1 in the fully turbulent flow region. In the 
streamwise direction, the grid density is high near the 
leading edge to resolve the large local flow gradients. 

Fig. (6) shows the results of the grid independent 
studies by continuously refining the grid in the stream- 
wise direction. The grids are 135 x 120, 215 x 120, 255 x 
120, 285 x 120, 335 x 120. There is no significant differ- 
ence in the skin friction coefficient, which is very sen- 
sitive to grid density, among those obtained from the 
255 x 120,285 x 120, and the 335 x 120 meshes. 

Fig.(7) shows the calculated skin friction coefficient 
using grids of 255 x 120 and 255 x 160. The results 
show little difference. These results indicate that the 
255 x 120 grid can provide a solution of a fairly high 
level of grid independency. 

• Model Comparisons : T3A 

The free stream turbulence intensity is 3% for case 
T3A. The computational mesh is 255 x 120. The inlet 
conditions and the mesh for the present model are iden- 
tical to that for the LS model. The Fig. (8) shows the 
comparison of computed and measured 1 skin friction 
coefficients, C/, as functions of Re x . The present model 
predicts very well for both the onset and the length of 
the transition region. The current model’s prediction of 
the variation of Cj through the transitional region also 
agree well with the measurement. The peak value of 
Cj, which occurs near the end of the transition zone, is 
also well predicted, indicating that the current model 


5 



can predict the location of the maximum momentum 
(or heat) transfer. 

The calculated and measured streamwise variations 
of the shape factor are shown in Fig. (9). The result 
from the LS model is also shown for comparison. 

• Model Comparisons : T3B 

For the T3B case, the free stream turbulence is 
about 6%. We have used a computational mesh of 
255 x 120 for both models. Again, the inlet conditions 
are the same for both the present model and the LS 
model. A comparison of the calculated and the mea- 
sured streamwise variation of Cj is shown in Fig. (10). 
The transition onset is predicted well by the present 
model. Compared with the data points available, the 
predicted length of the transition region appears to be 
somewhat shorter. The difference between the pre- 
dicted maximum value of Cj and the maximum value of 
C; of the available data is about 4%, which falls within 
the experimental margin of error for such a measure- 
ment. 

Fig. (11) shows the streamwise variation of the shape 
factors. The LS model predicts an earlier drop off than 
the measurement. The present model, on the other 
hand, predict a delayed occurrence of such a decrease 
of the shape factor. 

In Fig. (12), the calculated distributions of the dis- 
turbance kinetic energy, k, along the direction normal 
to plate at six streamwise locations are shown. The nor- 
mal distances y for the individual streamwise locations 
are nondimensionalized by their respective boundary 
layer thickness, S, and the k are nondimensionalized 
by the free stream velocity at the inlet. Note that the 
present model uses the local peak value of k to charac- 
terize the transition progress variable, /. At Re r = 10 4 , 
where the calculated flow is laminar, k varies almost lin- 
early within the boundary layer and there is no internal 
peak value. At later locations, an internal peak begins 
to appear in the k profile. The peak value increases 
as the flow goes through the transition process. At 
Rt z = 10 5 , the peak value of k appears to be leveling 
off and begin to decrease. According to the variation of 
peak value of fc, the model is predicting a gradually in- 
creasing turbulence intensity as the flow develop in the 
transition region. The turbulence intensity reaches a 
maximum level at about Re x = 10 5 , which roughly cor- 
responds to the calculated maximum C/. The results 
shows that the local peak value of k is an appropri- 


ate parameter for characterizing the bypass transition 
process. 

Summary 

A new model for bypass transitional flows has been 
described. The model is developed based on the obser- 
vation that the transitional boundary layers are high 
intermittent and, accordingly, an intermittency correc- 
tion function is incorporated in the turbulent eddy- 
viscosity. The model has been tested for bypass transi- 
tional flat plate boundary layers with different levels of 
turbulence intensity in the free stream using a Navier- 
Stokes equations solver. The results show that the 
present model can capture very well the effects of free 
stream turbulence on bypass transition. It should be 
noted that the Launder and Sharma’s turbulence model 
is not a dedicated bypass transition model. Therefore, 
its results shown here for the bypass transition cases 
should be viewed accordingly. 
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Fig.(7) Grid independent study with grid refinement 
in the direction normal to the plate. 



Fig. (8) Comparison of skin friction coefficients. T3A. 



Fig.(9) Comparison of shape factor distributions. T3A. 



Fig.(10) Comparison of skin friciton coefficients. T3B. 



Fig.(ll) Comparison of shape factor distributions.T3B. 



Fig. (12) Disturbance kinetic energy distrubutions in 
the direction normal to the fiat plate at different 
streamwise locations. T3B 
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